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, Abstract 

\Q • A preliminary group classification of the class 2D nonlinear heat equations ut = f(x,y,u,u x ,u y )(u xx + u yy ), where / is 
£SJ " arbitrary smooth function of the variables x, y, u, u x and u y using Lie method, is given. The paper is one of the few applications 
of an algebraic approach to the problem of group classification: the method of preliminary group classification. 
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1 Introduction 



It is well known that the symmetry group method plays an important role in the analysis of differential equations. 
The history of group classification methods goes back to Sophus Lie. The first paper on this subject is [1], where Lie 
» ^ proves that a linear two-dimensional second-order PDE may admit at most a three-parameter invariance group (apart 
\q \ from the trivial infinite parameter symmetry group, which is due to linearity). He computed the maximal invariance 
[ — . group of the one-dimensional heat conductivity equation and utilized this symmetry to construct its explicit solutions. 
' Saying it the modern way, he performed symmetry reduction of the heat equation. Nowadays symmetry reduction 
\ is one of the most powerful tools for solving nonlinear partial differential equations (PDEs). Recently, there have 
■ been several generalizations of the classical Lie group method for symmetry reductions. Ovsiannikov [2] developed 
' the method of partially invariant solutions. His approach is based on the concept of an equivalence group, which is a 
C^} , Lie transformation group acting in the extended space of independent variables, functions and their derivatives, and 
^ ■ preserving the class of partial differential equations under study. In an attempt to study nonlinear effects Saied and 
Hussain [3] gave some new similarity solutions of the (l+l)-nonlinear heat equation. Later Clarkson and Mansfield [4] 
studied classical and nonclassical symmetries of the (1+1) -heat equation and gave new reductions for the linear heat 
' equation and a catalogue of closed-form solutions for a special choice of the function f(x,y,u,u x ,u y ) that appears 
| in their model. In higher dimensions Servo [5] gave some conditional symmetries for a nonlinear heat equation while 
Goard et al. [6] studied the nonlinear heat equation in the degenerate case. Nonlinear heat equations in one or higher 
dimensions are also studied in literature by using both symmetry as well as other methods [7,8]. 
There are a number of papers to study (l+l)-nonlinear heat equations from the point of view of Lie symmetries 
method. The (2+l)-dimensional nonlinear heat equations 

U t = f(u)(U XX +Uyy), (1.1) 

are investigated in [9] and in present paper we studied 

U t = f(X, l/, U, U X ,Uy)(U XX + Uyy), (l .2) 

Similarity techniques are applied in [10,11,12,13] for (2+l)-dimensional wave equations. 
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2 Symmetry Methods 

Let a partial differential equation contains p dependent variables and q independent variables. The one-parameter 
Lie group of transformations 

Xi i — >Xi + eC (x,u) + 0(e 2 ); u a i — > u a + etp a (x, u) + 0(e 2 ), (2.3) 

where i = 1, . . . ,p and a = 1, . . . ,q. The action of the Lie group can be recovered from that of its associated 
infinitesimal generators, we consider general vector field 

i— 1 a— 1 

on the space of independent and dependent variables. The symmetry generator associated with (2.4) given by 

d d d d 

X = i 1 {x,y,t,u)—+S, 2 {x,y,t,u)—+i 3 {x,y,t,u)— + tp(x,y,t,u) — . (2.5) 

The second prolongation of X is the vector field 

X^> = X + (p— + <p u — + <p— + <p —— + l py—— + ip +(p yy^— + ip f +ip 



du x du y du t du xx du xy du xt du yy du yt du t t 

(2.6) 

that its coefficients are obtained with following formulas 

tp x = D x ip - u x D x £ - u y D x £ 2 - u t D x f, ip y = Dytp-UxDyZ 1 - u y D y ^ 2 -u t D y f, (2.7) 

<p* = D t cp - u x D t ^ - u v D t £? ~ u t B t f, tp xx = D x ip x - u xx D x £ - u xy D x f - u xt D x £ 

<p yy = Dy^ - U X yDy£} - UyyDy^ 2 ~ UyfDyf = D t </ ~ U^D^ ~ Uy t D t ^ ~ U„ Df^ 

(fi Vt = Dyip 1 ~ U X yDy£} - UyyDyi 2 ~ U yt D y ^ if^ = D ^ ~ U xt D t ^ ~ U^D^ ~ U tt 

where the operators D x , D y and D t denote the total derivatives with respect to x,y and t: 
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d d d d d 

— + U x — ~t~ I" ^xy~^. ~t~ V>xt~Z. 

ox ou ou x ou y ou t 

d d d d d 



Dy= dy- +U ^ + U ™du~y + U ^du~ x +U ^du- t + --- (2 ' 8) 

d d d d d 

D t = — + u t — + u tt - h u tx - h u ty - \- ... 

at ou ou t ou x Ouy 

By theorem 6.5. in [14], X^[ut — f(x,y,u,u x ,u y )(u xx + u yy )]\^ — whenever 

U t ~ f(X, V, U, U X , Uy)(U XX + Uyy) = 0. (2.9) 

Since 

X ( - 2) [u t -f(x,y,u,u x ,Uy)(u xx +Uyy)} = Lp t -(f x ^ 1 +f y ^ 2 +f u (p+f Ux (f x +f Uy Lp y )(u xx +u yy )-f(x,y,u,u x ,u y )((f xx 
therefore we obtain the following determining function: 

V* ~ (fxt 1 + fyt 2 + fu<P + fu x V x + fu y f y )(u xx + u yy ) - f(x,y,u,u x ,u y )(ip xx + ip^) = 0. (2.10) 
In the case of arbitrary / it follows 

e=e = v=o, (2.H) 
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or 

e = e = ^ = o, e = c. (2.12) 

Therefore, for arbitrary f(x,y,u,u x ,u y ) Eq. (1.1) admits the one-dimensional Lie algebra gi, with the basis 

* = I (2.13) 

Qi is called the principle Lie algebra for Eq. (1.1). So, the remaining part of the group classification is to specify the 
coefficient / such that Eq. (1.1) admits an extension of the principal algebra gi. Usually, the group classification is 
obtained by inspecting the determining equation. But in our case the complete solution of the determining equation 
(2.10) is a wasteful venture. Therefore, we don't solve the determining equation but, instead we obtain a partial group 
classification of Eq. (1.1) via the so-called method of preliminary group classification. This method was suggested in 
[10] and applied when an equivalence group is generated by a finite-dimensional Lie algebra Qg. The essential part 
of the method is the classification of all nonsimilar subalgebras of gg. Actually, the application of the method is 
simple and effective when the classification is based on finite-dimensional equivalence algebra Qg. 

3 Equivalence transformations 

An equivalence transformation is a nondegenerate change of the variables t,x,y,u taking any equation of the form 
(1.1) into an equation of the same form, generally speaking, with different f(x, y, u, u x , u y ). The set of all equivalence 
transformations forms an equivalence group $ . We shall find a continuous subgroup $c of it making use of the 
infinitesimal method. 

We consider an operator of the group $c m the form 

d d d d d 

y = t 1 (x,y,t,u)— + £, 2 (x,y,t,u)—+ £ 3 (x,y,t,u)— + <p(x,y,t,u)— + fj,(x,y,t,u,u x ,u y ,u t , f) — , (3.14) 

from the invariance conditions of Eq. (1.1) written as the system: 

ut ~ f(x,y,u,u x ,u y )(u xx +u yy ) = 0, (3.15) 
ft = fu t = 0, 

where u and / arc considered as differential variables: u on the space (x,y,t) and / on the extended space 
(x,y,t,u,u x ,u y ). 

The invariance conditions of the system (3.15) are 

Y {2) (u t - f(x,y,u,u x ,u v )(u xx +u yy )) = 0, (3.16) 
y(2) (/t) = y (2) (/ut) = 0; 

where is the prolongation of the operator (3.14): 

au x du y au t ou xx du xy du xt ou yy 

+ v vt ^- + f tt ^- + + m"* — ~~ — 



diiyt du tt df t df Ut ' 

The coefficients ip x , (p v , (p*, (p xx , tp xy , ip xt , (p vv , ip yt , (p u are given in (2.7) and the other coefficients of (3.17) are obtained 
by applying the prolongation procedure to differential variables / with independent variables (x,y,t,u,u x ,u y ,ut). 
we have 

fjL* = D t in) - f x D t {?) - fyDtif) - fMv) ~ fuM^J ~ fu y D t (v v \ (3-18) 
M"' =DM - f x D Ut (e) - fy£>u t (f) - fuD Ut (<p) - fu x D Ut ( v x ) - f Uy D Ut (<pV), (3.19) 
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where 



So, we have the following prolongation formulas: 

- IH - Ul - fy& - fuVt - In, {<f?)t - fu y {<P V )U (3.21) 

By the invariancc conditions (3.16) give rise to 

H* = n u * = 0, (3.22) 
that is hold for every /. Substituting (3.22) into (3.21), we obtain 

\H = A*« t = 

e x =tf=<Pt = (3.23) 

{<f?) t = (<p x ) ut = {ip») t = (<py) ut = o 

Moreover with substituting (3.17) into (3.16) we obtain 

ip* - f(x, y, u, u x , u y )((p xx + ip yy ) - fi(u xx + u yy ) = 0. (3.24) 

We are left with a polynomial equation involving the various derivatives of u(x,y,t) whose coefficients are certain 
derivatives of £ 2 , £ 3 and <p. Since £ 2 , £ 3 , <^ only depend on x, y, t, u we can equate the individual coefficients to 
zero, leading to the complete set of determining equations: 

(3.25) 
(3.26) 
(3.27) 
(3.28) 
(3.29) 
(3.30) 
(3.31) 
(3.32) 
(3.33) 

so, we find that 

^(x) = Cl x + c 2 y + c 3 , i 2 {t)=c iy + Ci , £ 3 (t)=a(t), 

<p(x,y,u) = au + 0(x,y), n = (ci - a'(t))f, (3.34) 



e 




e 


= e(y) 


e 




Vuu 


= 


2<Pxu 


^xx 1 ^yy 


<Pyu 




Wu 




A* 






= f(<Pxx + <Pyy) 



with constants c\, c 2 , c 3 and C4, also we have /3 XX = — /3, 



;/?/• 



We summarize: The class of Eq. (1.2) has an infinite continuous group of equivalence transformations generated 
by the following infinitesimal operators: 

d d d d d 

Y = ( Cl x + c 2 y + c 3 )— + [ay + c 4 )— + a(t)— + {cm + 0(x, y))-^ + (ci - a'{t))f—. (3.35) 
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Table 1 

Commutation relations satisfied by infinitesimal generators in (4.38) 



[,] 


Yi 


Y 2 


Y 3 


Yi 


Y 5 


Y 6 


Yi 








-Y 3 


-Y 4 





-Y B 


Y 2 











-Y 3 








Y 3 


Y 3 

















Y 4 


Y 4 


Y 3 














Y 5 




















Y 6 


Y 6 


















Table 2 

Adjoint relations satisfied by infinitesimal generators in (4.38) 



[,] 


Yi 


Y 2 


Y 3 


Y 4 


Y 5 


Y 6 


Yi 


Yi 


Y 2 


e s Y 3 


e a Y 4 


Y 5 


e a Y 6 


Y 2 


Yi 


Y 2 


Y 3 


Y 4 + sY 3 


Y 5 


Y 6 


Y 3 


Yi - sY 3 


Y 2 


Y 3 


Y 4 


Y 5 


Y 6 


Y 4 


Yi - sY 4 


Y 2 - sY 3 


Y 3 


Y 4 


Y 5 


Y 6 


Y 5 


Yi 


Y 2 


Y 3 


Y 4 


Y 5 


Y 6 


Y 6 


Yi - sY 6 


Y 2 


Y 3 


Y 4 


Y 5 


Y 6 



Therefore the symmetry algebra of the Burgers' equation (1.2) is spanned by the vector fields 

d d d d d x „ 5 T „ 3 ^ r d . 

Yi =x— +y— +t-+u— + f— , Y 2 =y—, Y z = —, Y 4 = — 3.36 

ox oy ot ou of ox ox oy 

Moreover, in the group of equivalence transformations there are included also discrete transformations, i.e., reflections 
t — > -t, x — > -x, u — > -u, f — ► -/. (3.37) 



4 Preliminary group classification 

One can observe in many applications of group analysis that most of extensions of the principal Lie algebra admitted 
by the equation under consideration are taken from the equivalence algebra Qg. We call these extensions S -extensions 
of the principal Lie algebra. The classification of all nonequivalent equations (with respect to a given equivalence 
group Gg,) admitting ^-extensions of the principal Lie algebra is called a preliminary group classification. Here, 
Gg is not necessarily the largest equivalence group but, it can be any subgroup of the group of all equivalence 
transformations. 

So, we can take any finite-dimensional subalgebra (desirable as large as possible) of an infinite-dimensional algebra 
with basis (3.30) and use it for a preliminary group classification. We select the subalgebra g e spanned on the 
following operators: 

d d d d , d (9 d () 

Yl= X 7T + V-^- +t ~^l +u 7r + F 2=y^-, Y 3 = —, Y 4 = t~, 

Ox oy ot ou of ox ox oy 

*=l-4' 14 <4 - 38 > 
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The communication relations between these vector fields is given in Table 1. To each s-paramctcr subgroup there 
corresponds a family of group invariant solutions. So, in general, it is quite impossible to determine all possible 
group-invariant solutions of a PDE. In order to minimize this search, it is useful to construct the optimal system of 
solutions. It is well known that the problem of the construction of the optimal system of solutions is equivalent to 
that of the construction of the optimal system of subalgebras [2,12]. Here, we will deal with the construction of the 
optimal system of subalgebras of £5. 

Let G be a Lie group, with q its Lie algebra. Each element TeG yields inner automorphism T a — ► TT a T~ x of the 
group G. Every automorphism of the group G induces an automorphism of 0. The set of all these automorphism is 
a Lie group called the adjoint group G A . The Lie algebra of G A is the adjoint algebra of 0, defined as follows. Let 
us have two infinitesimal generators X, Y G L. The linear mapping AdX(Y) : Y — ► [-^jY] is an automorphism of 
0, called the inner derivation of g. The set of all inner derivations a,dX(Y)(X, Y E g) together with the Lie bracket 
[AdX, AdF] = AdLY, Y] is a Lie algebra g A called the adjoint algebra of g. Clearly g A is the Lie algebra of G A . Two 
subalgebras in g are conjugate (or similar) if there is a transformation of G A which takes one subalgcbra into the 
other. The collection of pairwise non-conjugate s-dimensional subalgebras is the optimal system of subalgebras of 
order s. The construction of the one-dimensional optimal system of subalgebras can be carried out by using a global 
matrix of the adjoint transformations as suggested by Ovsiannikov [2]. The latter problem, tends to determine a list 
(that is called an optimal system) of conjugacy inequivalent subalgebras with the property that any other subalgebra 
is equivalent to a unique member of the list under some element of the adjoint representation i.e. f) Ad(g) f) for some 
g of a considered Lie group. Thus we will deal with the construction of the optimal system of subalgebras of g 6 - 

The adjoint action is given by the Lie series 

Ad(exp( S Yi))Yj =Y 3 -s [Y t ,Y 3 ] + *- [Y u [Y h Y 3 ]] , (4.39) 

where s is a parameter and i,j = 1, • • • ,6. The adjoint representations of 06 is listed in Table 2, it consists the 
separate adjoint actions of each element of 6 on all other elements. 

Theorem 4.1. An optimal system of one- dimensional Lie subalgebras of general Burgers' equation (1-2) is provided 
by those generated by 



1) Y 1 -- 


= Yi = xd x + yd y + td t + ud u + fd f , 


2) y 2 = 


Y 2 = yd x , 








3) Y 3 - 


= ~Yi = -By, 


4) y 4 = 


Y!+Y 5 = 


xd x 4 


-yd y + (t+l)d t 


+ ud u , 


5) Y 5 -- 


= Yi - y 2 = (x - y)d x + ydy + td t + ud u - 


\-fdf, 6)y 6 = 


Y 2 -Y 4 = 


ydx - 


dy, 




7) Y 7 -- 


= -y 4 + y 6 = -dy + d u , 


8) y 8 = 


-Yi - Y 6 : 


- -d - d 

— Vy U U , 




9) Y 9 -- 


= Y 2 + Y 5 = yd x + d t ~ fd f , 


10) y 10 = 


y 2 - y 5 = 


■ yd x - 


- dt + fdf, 




11) Y 11 


= Y 2 + Y 6 = yd x + d u , 


12) y 12 = 


y 2 - y 6 - 


yd x - 


- d u , 




13) y 13 


= Yi+Y 2 = (x + y)d x + ydy + td t + ud u 


+ fdf, 14) y 14 = 


-Yi + Y 5 


+ Y 6 -- 


= -d y + d t + d u 


- fdf, 


15) y 15 


= Y 2 -Y 4 -Y 5 + Y 6 = yd x - d y - d t + 9, 


> + fdf, 16) y 16 - 


Y 2 -Y 4 + 


Y 6 = 


yd x -d y + d u , 




17) y 17 


= Y 2 -Y 4 + Y 5 -Y 6 = yd x - d y + d t - 9, 


fdf, is) y 18 = 


Y 2 -Y 4 - 


Y 6 = 


yd x -d y - d u , 




19) y 19 


= Y 1 +Y 2 + Y 5 = (x + y)d x + (t+ l)d t + ud u , 20) Y 2a = 


-Y 2 + Y 5 + Y 6 = 


yd x + d t + d u - 


fdf, 


21) y 21 


= Y 2 + Y 5 -Y 6 = yd x + d t -d u - fd f , 


22) y 22 = 


Y 2 -Y 5 - 


Y 6 = 


yd x -d t -d u + 


fdf, 


23) y 23 


= Y 2 -Y 5 + Y 6 = yd x -d t + d u + fdf, 


24) y 24 - 


^ -Yi - y 5 


-Y 6 


= -d y - d t - d u 


+ fdf, 


25) y 25 


= -Yi - y 5 + y 6 = -dy -d t + d u + fd f , 


26) y 26 = 


-Yi + y 5 


-Y e 


= -dy + d t -d u 


- fdf, 


27) y 27 


= y 2 -y 4 + y 5 + y 6 = y d x -d y + d t + a, 












28) y 28 


= Y 1 +Y 2 -Y 5 = {x + y)d x + ydy + (t - 


i)d t + udu + 2fd f 










29) y 29 


= Y 1 -Y 2 -Y 5 = {x- y)d x + ydy + (t - 


l)d t + ud u + 2/9/ 










30) y 31 


= Y 1 -Y 2 + Y 6 = (x- y)d x + yd y + (t + 


l)9 t + ud u , 










31) y 31 


= Yi-Y 5 = xd x + ydy + (t - l)d t + ud u 


+ Vdf 










32) y 32 


= Y 2 - Y 4 - y 5 - y 6 = yd x - dy- d t - 9, 


, + fdf, 
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Proof. Let Qe is the symmetry algebra of Eq. (1.2) with adjoint representation determined in Table 2 and 

Y = aiFi + a 2 Y 2 + a 3 Y 3 + a 4 Y 4 + a 5 Y 5 + a 6 Y 6 , (4.40) 

is a nonzero vector field of g@. We will simplify as many of the coefficients c^; i = 1, . . . , 6, as possible through proper 
adjoint applications on Y. We follow our aim in the below easy cases: 
Case 1: 

At first, assume that a\ 7^ 0. Scaling Y if necessary, we can assume that a\ = 1 and so we get 

Y = Yi + a 2 Y 2 + a 3 Y 3 + a 4 Y 4 + a 5 Y 5 + a 6 Y 6 . (4.41) 
Using the table of adjoint (Table 2) , if we act on Y with Ad(exp(a3Y3)), the coefficient of Y 3 can be vanished: 

Y' = Yi + a 2 Y 2 + a 4 Y 4 + a 5 Y 5 + a 6 Y 6 . (4.42) 
Then we apply Ad(exp(a4Y4)) on Y' to cancel the coefficient of Y 4 : 

Y" = Y 4 + a 2 Y 2 + a 6 Y 6 + a 6 Y e . (4.43) 
At last, we apply Ad(exp(a6Yes)) on Y" to cancel the coefficient of Y e : 

Y'" = Y x + a 2 Y 2 + a 5 Y 5 . (4.44) 
Case la: 

If a 2 , CI5 7^ then we can make the coefficient of Y 2 and Y 5 either +1 or —1. Thus any one-dimensional subalgebra 
generated by Y with a 2 , a$ ^ is equivalent to one generated by Y\ ± Y 2 ± Y 5 which introduce parts 19), 28), 29) 
and 30) of the theorem. 
Case lb: 

For a 2 — 0, a 5 7^ we can see that each one-dimensional subalgebra generated by Y is equivalent to one generated 
by Yi ± Y 5 which introduce parts 4) and 31) of the theorem. 
Case lc: 

For a 2 7^ 0, a 5 = 0, each one-dimensional subalgebra generated by Y is equivalent to one generated by Y\ ± Y 2 which 
introduce parts 5) and 13) of the theorem. 
Case Id: 

For a 2 — 0, a$ = 0, each one-dimensional subalgebra generated by Y is equivalent to one generated by Y\ which 
introduce parts 1) of the theorem. 

Case 2: 

The remaining one-dimensional subalgebras are spanned by vector fields of the form Y with ai = 0. 
Case 2a: 

If a 4 7^ then by scaling Y, we can assume that a 4 = —1. Now by the action of Ad(expa3Y3)) on Y, we can cancel 
the coefficient of Y 3 : 

Y = a 2 Y 2 - Y 4 + a 5 Y 5 + a e Y 6 . (4.45) 
Let a 2 7^ then by scaling Y, we can assume that a 2 = 1, and we have 

Y' =Y 2 -Y 4 + a 5 Y 5 + a 6 Y e . (4.46) 
Case 2a- 1: 

Suppose d5 = d6 = 0, then the one-dimensional subalgebra generated by Y is equivalent to one generated by Y 2 — Y 4 
which introduce parts 6). 
Case 2a-2: 

Suppose a 5 = 0, a 6 7^ 0, all of the one-dimensional subalgebra generated by Y is equivalent to one generated by 
Y 2 — Y 4 ± Ye; which introduce parts 16) and 18). 
Case 2a-3: 

Suppose a 5 7^ 0, a 6 7^ 0, all of the one-dimensional subalgebra generated by Y is equivalent to one generated by 
Y 2 -Y 4 ±Y 5 ± Y 6 which introduce parts 15), 17), 27), and 32). 
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Now if a 2 = 0, we have 



Y" = -Y i + a 5 Y 5 + a 6 Y e . (4.47) 
Case 2a- 4: 

Suppose a 5 = a 6 = 0, then the one-dimensional subalgebra generated by Y is equivalent to one generated by —Y4 
which introduce parts 3). 
Case 2a-5: 

Suppose as = 0, ae ^ 0, all of the one-dimensional subalgebra generated by Y is equivalent to one generated by 
—Y4 ± Y 6 which introduce parts 7) and 8). 
Case 2a-6: 

Suppose a 5 ^ 0, a 6 7^ 0, all of the one-dimensional subalgebra generated by Y is equivalent to one generated by 
-Y A ±Y 5 ±Y 6 which introduce parts 14), 24), 25) and 26). 

Case 2b: 

Let 04 = then Y is in the form 

Y = a 2 Y 2 + a 5 Y 5 + a 6 Y 6 . (4.48) 
Suppose that a 2 7^ then if necessary we can let it equal to 1 and we obtain 

?' = Y 2 + a 5 Y 5 + a 6 Y 6 . (4.49) 
Case 2b- 1: 

Let (Z5 = a% = 0, then Y 2 is remained and find 2) section of the theorem. 
Case 2b-2: 

If a 5 7^ 0, a 6 7^ 0, then Y' is equal to Y 2 ± Y 5 ± Y 6 - Hence this case suggests part 20), 21), 22) and 23). 
Case 2b-3: 

If 05 7^ 0, aa — 0, then Y' = Y 2 ±Y 5 . Hence this case suggests part 9) and 10). 
Case 2b-4: 

If as = 0, a 6 7^ 0, then Y 2 ± Y & is obtained. So, this case suggests part 11) and 12). 

There is not any more possible case for studying and the proof is complete. □ 

The coefficients / of Eq. (1.2) depend on the variables x,y,u,u X7 u y . Therefore, we take their optimal system's 
projections on the space (x, y, u, u x , u y , /). we have 

1) Z 1 = Y 1 = xd x + yd y + ud u + fd f , 17) Z 17 = Y 17 = (x - y)d x + yd y + ud u + 2/9/, 

2) Z 2 = Y 2 = yd x , 18) Z 1S - F 18 = {x + y)d x + ud Ul 

3) Z 3 = Y 3 = -d y , 19) Z 19 = Y 19 = yd x + d u - fd f , 

4) Z 4 = F 4 = (x + y)d x + yd v + ud u + fd f , 20) Z 2a = Y 20 = yd x - d u - fd f , 

5) Z 5 = Y 5 = xd x + ydy + ud u , 21) Z 21 = Y 21 = yd x - d u + fdf, 

6) Z 6 = Y 6 = xd x + ydy + ud u + 2fd f , 22) Z 22 = Y 22 = yd x + d u + fd f , 

(4.50) 

7) Z 7 = Y 7 = (x - y)d x + ydy + ud u + fd f , 23) Z 23 = Y 23 - yd x - d y + d u , 

8) Z S =Y S = yd x - d y , 24) Z 24 = Y 24 = yd x - d v - d u , 

9) z 9 = y 9 = -d y + d u , 25) z 25 = r 25 = -d y + d u - fd f , 

10 z 10 = y w = -dy - d u , 26) z 26 = r 26 = -dy -d u + fd f , 

11 Z 11 = Y 11 = yd x - fdf, 27) Z 27 = Y 27 = -d y + d u + fdf, 
12) Z 12 = Y 12 = yd x + fdf, 28) Z 28 = Y 28 = -d y - d u - fdf, 
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13) Z 13 


_ yl3 


= + 9 U , 


29) Z 29 


_ y29 




-8y 


+ d u 


~fdf, 


14) Z 14 


_ yl4 


= yd x - d u , 


30) Z 30 


_ y30 


= yd* 


-dy 


- d u 


+ fdf, 


15) Z 15 


= y 15 


= (x - y)d x + yd v + ud u , 


31) Z 31 


_ y31 


= y<9 x 


- dy 


+ d u 


+ fdf, 


16) Z 16 


_ y!6 


= {x + y)d x + yd y + ud u + 2fd f , 


32) Z 32 


_ y32 


= yd x 


-dy 


- d u 


-fdf, 



(4.51) 



Proposition 4.2. LetQ m := (Yi, . . . ,Y m ), be an m- dimensional algebra. Denote byY l (i = 1, . . . , r, < r < m, r G N) 
an optimal system of one- dimensional subalgebras of Q m and by Z l (i = 1, • • • ,(,0<(<r,teN) £/ie projections of 
Y l , i.e., Z l — pr(F J ). // equations 

f = $(x,y,u,Ux,u y ), (4.52) 
are invariant with respect to the optimal system Z l then the equation 

u t = y, u, u x ,u y )(u xx + u yy ), (4.53) 
admits the operators X 1 = projection of Y l on (t,x,y,u,u x ,u y ). 
Proposition 4.3. Let Eq. (4-53) and the equation 

U t = $'(X, y, U, U X , Uy)(U XX + Uyy), (4.54) 

be constructed according to Proposition 4-2. via optimal systems Z l and Z 1 ' , respectively. If the subalgebras spanned 
on the optimal systems Z l and Z % , respectively, are similar in Q m , then Eqs. (4-53) and (4-54) are equivalent with 
respect to the equivalence group G m , generated by Q m . 

Now we apply Proposition 4.2. and Proposition 4.3. to the optimal system (3.16) and obtain all nonequivalcnt Eq. 
(1.2) admitting ^-extensions of the principal Lie algebra g, by one dimension, i.e., equations of the form (1.2) such 
that they admit, together with the one basic operators (4.55) of g, also a second operator X^ 2 \ For every case, when 
this extension occurs, we indicate the corresponding coefficients /, g and the additional operator X^ 2 \ 

We perform the algorithm passing from operators Z l (i = 1, • • • , 32) to / and X^ via the following example. 
Let consider the vector field 

Z 32 =yd x -d y -d u - fd f , (4.55) 

then the characteristic equation corresponding to Z 6 is 

dx dy du df , , „. 

7 = ^ = n = z? < 4 - 56 > 

and can be taken in the form 

h=u+-, h = e%f. (4.57) 

y 

From the invariance equations we can write 

h = (4.58) 
it follows that 

/ = e -f$(A), (4.59) 
where A = I\. 
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From Proposition 4.2. applied to the operator Z & we obtain the additional operator 

yd x - d y + 8 t - d u . (4.60) 

After similar calculations applied to all operators (3.16) we obtain the following result (Table 3) for the preliminary 
group classification of Eq. (1.2) admitting an extension g 3 of the principal Lie algebra jji. 

5 Conclusion 

In this paper, following the classical Lie method, the preliminary group classification for the class of heat equation 
(1.2) and investigated the algebraic structure of the symmetry groups for this equation, is obtained. The classification 
is obtained by constructing an optimal system with the aid of Propositions 4.2. and 4.3.. The result of the work is 
summarized in Table 3. Of course it is also possible to obtain the corresponding reduced equations for all the cases 
in the classification reported in Table 3. 
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Table 3 

The result of the classification 



N 


Z 


Invariant 


Equation 


Additional operator 


1 


Z 1 


u 

X 


u t 


= X<&(U XX + Uyy) 


xdx + yd y + td t + ud u 


2 


Z 2 


U 


u t 


= ®{UxX + Uyy) 


ydx 


3 


Z 3 


U 


Ut 


= &{U XX + Uyy) 


-8y 


4 


Z 4 


u 


Ut 


= y<b{U X X + Uyy) 


(x + y)d x + ydy + ud u 


5 


z 5 


u 


Ut 


= ${UxX + Uyy) 


xdx + ydy + (t + l)d t + ud u 


6 


z 6 


u 


Ut 


= X 2 $(U X X + Uyy) 


xdx + ydy + (t - l)d t + udu 


7 


z 7 


u 


Ut 


= (X - y)${U X X + Uyy) 


(x - y)dx + ydy + td t + ud u 


8 


z s 


U 


Ut 


= ®{UxX + Uyy) 


ydx - d y 


9 


z 9 


U 


Ut 


= ®{UxX + Uyy) 


-d y + d u 


10 




X 


Ut 


= Q(UxX + Uyy) 


-dy - d u 


11 


z 11 


U 


Ut 


= e + Uyy) 


ydx + d t 


12 


Z 12 


u 


Ut 


= e»$(tt« + Uyy) 


ydx - d t 


13 


Z 13 


u- | 


Ut 


= ®{UxX + Uyy) 


ydx + du 


14 


Z 14 


u+ | 


Ut 


= ®{UxX + Uyy) 


ydx - du 


15 


Z 18 


u 


Ut 


= Q(UxX + Uyy) 


(x - y)d x + ydy + (t + l)d t + ud u 


16 


Z 16 


u 


Ut 


= (X + y) 2 <&{U XX + Uyy) 


(x + y)dx + ydy + (t - l)d t + ud u 


17 


z 17 


u 


Ut 


= (X -y) 2 &(U XX +Uyy) 


(x - y)d x + ydy + (t - l)d t + ud u 


18 


Z 18 


u 


Ut 


= Q(UxX + Uyy) 


(x + y)dx + {t+ l)d t + udu 


19 


Z 19 


u ~-y 


u t 


= e~y$(u XX + Uyy) 


ydx + d t + d u 


20 


Z 20 


u+l 


Ut 


= e _ »$(u iX + Uyy) 


ydx + d t - du 


21 


Z 21 


u+ | 


Ut 


= + Uyy) 


ydx - d t - d u 


22 


Z 22 


u- | 


Ut 


= e y$(U X X + Uyy) 


ydx - d t + d u 


23 


Z 23 


u- | 


Ut 


= ®{UxX + Uyy) 


ydx - d y + d u 


24 


Z 24 


u+ f 


Ut 


= e V $(Uxx + Uyy) 


ydx - d y - d u 


25 


Z 25 


M + y 


Ut 


= ${U X X + Uyy) 


-d y + d t + d u 


26 


Z 26 


u-y 


Ut 


= e~ V $(Uxx + Uyy) 


-d y - d t - d u 


27 


Z 27 


u + y 


Ut 


= e~ V $(UxX + Uyy) 


-d y - d t + du 


28 


Z 28 


u-y 


Ut 


= e V $(Uxx + Uyy) 


-d y + d t - d u 


29 


Z 29 


y 


Ut 


= $(Uxx + Uyy) 


ydx - dy + dt + du 


30 


Z 30 


y 


Ut 


= e»$(ll„ + Uyy) 


ydx - d y - d t - d u 


31 


Z 31 


u- £ 
a 


Ut 


= eV$(UxX + Uyy) 


ydx - d y - d t + du 


32 


Z 32 


u+ - 

y 


Ut 


= e~y$(Uxx + Uyy) 


ydx - d y + d t - d u 
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